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Abstract —Pilot contamination is known to be one of the main 
impairments for massive MIMO multi-cell communications. In¬ 
spired by the concept of fractional frequency reuse and by recent 
contributions on pilot reutilization among non-adjacent cells, we 
propose a new pilot allocation scheme to mitigate this effect. The 
key idea is to allow users in neighboring cells that are closest to 
their base stations to reuse the same pilot sequences. Focusing 
on the uplink, we obtain expressions for the overall spectral 
efficiency per cell for different linear combining techniques at 
the base station and use them to obtain both the optimal pilot 
reuse parameters and the optimal number of scheduled users. 
Numerical results show a remarkable improvement in terms of 
spectral efficiency with respect to the existing techniques. 

Index Terms —Channel estimation, fractional frequency reuse, 
fractional pilot reuse, massive MIMO, pilot contamination. 

I. Introduction 

Massive multiple-input multiple-output (MIMO) refers to 
a type of cellular network based on multiuser MIMO in 
which the number of antennas at the base station (BS) is 
much larger than the number of served user terminals Q-©. 
Under certain propagation conditions, such a system renders 
quasi-orthogonal channels among users and very simple linear 
processing is shown to be optimal in terms of throughput. 

In this context, a problem that appears when obtaining 
channel state information is that the number of available pilot 
sequences is finite and limited by the channel behavior, as 
their duration cannot span larger than the coherence interval 
of the channel to be estimated. In consequence, mobile users 
in different cells might have to reuse the same pilot sequences, 
thus resulting into corrupted channel estimates at each BS: this 
is known as the pilot contamination effect © and represents a 
major impairment affecting massive MIMO communications. 

This problem has been intensively investigated in the litera¬ 
ture in the last few years. Eor instance, © proposed a low-rate 
coordination phase between different BSs; © introduced an 
elaborated scheme for OEDM-based massive MIMO systems, 
with differentiated downlink and uplink scheduled training 
phases; © proposed a subspace projection to improve the 
channel estimation accuracy; and © also exploited coordi¬ 
nation, comparing the amount of training needed with respect 
to the uncoordinated case and its implications. 

More recently, the idea of pilot reuse was introduced in ©. 
In this respect, |T§, |TT) advocated for employing orthogonal 
pilot subsets in adjacent cells, optimizing the number of 
required subsets and the number of scheduled users per cell 
that maximizes the overall spectral efficiency, for both uplink 
and downlink. This idea is borrowed from the well-established 


concept of frequency reuse in cellular systems p2) and is 
therein coined as “fractional pilot reuse”: given an integer 
/3 > 1 representing the number of orthogonal pilot subsets 
(or, equivalently, the ratio between the number of available 
sequences and the number of scheduled users per cell), for 
each cell, only a fraction 1//3 of the interfering cells reuse the 
same pilot subset. 

On the other hand, the concept of fractional frequency 
reuse (EER) consists in splitting the available bandwidth into 
a + 1 frequency sub-bands, where the integer a > 1 specifies 
the frequency reuse factor among the cell edges and the 
additional sub-band is reused in all the cell centers (see m 
and references therein). The difference between traditional 
frequency reuse and EER is illustrated in Eigure [T] 

The idea behind this paper is to adapt the aforementioned 
concept of EER to the pilot domain and assign the pilot 
sequences to the users within the cells in the same way as 
different frequency sub-bands are allocated in EER. Under 
these premises, we propose a novel channel estimation scheme 
where a fraction of users within each cell reuses the same pilot 
subset across the whole system, while the rest are allocated 
orthogonal subsets depending on a reuse parameter. We also 
refer to this approach as fractional pilot reuse, but in the 
sense that a fraction of the users per cell reuse the same 
pilot subset in all the cells. In this regard, we derive closed- 
form expressions for the spectral efficiency and analyze the 
performance in the large-antenna regime. Numerical results 
show that such deployment outperforms the proposal in | [T0| 
and offers remarkably good throughput, especially for large 
numbers of antennas at the BS. 

The rest of the paper is structured as follows. In Section [11] 
we introduce the system model on which we build our novel 
framework based on fractional pilot reuse in Section III In 
Section EYl numerical results are reported to corroborate the 
proposed scheme. Einally, Section |V] draws some concluding 
remarks. 


II. System Model 

Let us consider a massive MIMO cellular network over a set 
of cells denoted by C. The BS of each cell I G C is equipped 
with N antennas and serves K single-antenna users in the up¬ 
link; in this context, N is fixed while K is chosen adaptively. 
We model the channels as block-fading, where hj{zik) G 
represents the channel response between BS j and user k in 
cell I, with random coordinates G during its coherence 
block consisting of Tik channel uses: in the following, we 
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Figure 1. Traditional frequency reuse with a = 3 sub-bands and the 
equivalent FFR setup. 


consider that the coherence block of the whole system spans 
T = min/channel uses0 The channel realizations 
are drawn from a circularly symmetric complex Gaussian 
distribution, i.e., hj(z//j) ~ CAf(0, dj(z/fe)I/v), where dj{zii^) 
represents the variance of the channel attenuation from BS j 
to user k in cell I (cf. pT)): this term accounts for propagation 
impairments such as pathloss and shadowing. 

Let us introduce ^ x/^. = as the vec¬ 

tor signal transmitted by user k in cell I, normalized as 
E{|a:/fc(f)p} = 1, Vf S [IjT]. The received signal at BS j in 
a coherence block, denoted by Yj € is thus given by 

K 

Y, = E E (z/fc)x« + N, (1) 

l^C fc=l 


where pik is the transmit power allocated to user k in cell I and 
Nj € is the additive noise at BS j, whose elements are 

distributed independently as CAf(0,(T^). As in | [TI) , we adopt 
a statistic-aware power control^ where the signals from user k 
in cell I are allocated the power pik = p/di{zik), with p > 0. 

We take into account a pilot book V of B G [1,T] 
orthogonal pilot sequences dehned as 








B, if a = b 
0, otherwise 


( 2 ) 


with each V/ G C®; in the following, we will use iik to 
denote the index of the pilot sequence assigned to user k in 
cell 1. Therefore, the vector signal of user k in cell I can be 
written as x/^ = (u/fc, where U/^. G represents 

the information vector. We can thus decompose the received 
signal in ([T]) as Y, = (Y, Y,) with 


K 

y. = EE y^hj{zik)uff. + N 

IGC k^l 
K 

IGC /c=l 


+ N, 

(3) 

+ E 

(4) 


with Nj = (Nj Nj). Moreover, we assume that BS j applies 


user with high mobility, which would consistently decrease the overall 
coherence block T, can be seiwed in a different time-slot in order not to 
compromise the channel estimation of the other users. 

^Such power-control policy makes the average effective channel gain equal 
for all users. 



Figure 2. Example of fractional pilot reuse (FPR) with K = 5, 0 = 3, and 
= 2/K: in this setting, B = 11 orthogonal pilot sequences are required. 


a linear receive combining vector gjk G to the received 
signal, i.e., g^^Yj , so as to amplify the signal coming from its 
user k while rejecting the interference caused by the others. 
In this respect, we consider; 

a) maximum ratio combining (MRC), i.e., a traditional passive 
rejection scheme; 

b) pilot-based zero-forcing combining (P-ZFC), recently pro¬ 
posed in which actively suppresses both intra- and 
inter-cell interference. 

III. Fractional Pilot Reuse 
A. Problem Formulation 

Inspired by the concept of FFR, we split the pilot book V 
into /3 -f 1 subsets, where the integer /3 > 1 is referred to as 
integer reuse factor. The users that, inside each cell I, exhibit 
a low variance of the channel attenuation with respect to their 
BS (i.e., di{-)) reuse the /3 subsets in non-adjacent cells; on 
the other hand, the users that, inside each cell I, enjoy a high 
dif) are allocated the same (/3-|- l)-th subset across the whole 
system|^Hence, we dehne the new parameter /3f G [0,1) as the 
fractional reuse factor, which represents the fraction of users 
per cell that are assigned the additional pilot subset. We refer 
to this approach as fractional pilot reuse (FPR); an illustrative 
example is shown in Figure The number of pilot sequences 
(and, hence, their length) is determined as 

+ (l-/3f)/3). (5) 

Note that /3f = 1 would imply assigning the same B = K 
sequences to each cell with no pilot reuse whatsoever (this 
case is prevented by the dehnition of /3f). On the other hand, 

^In Section [T^ we consider a standard pathloss model in which the users 
are divided into two groups simply based on the distance from their BS. 
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if /3f = 0, the formulation trivially collapses to B = PK as 
in the baseline GD- 

Our objective is to find the pilot reuse parameters /3 and 
Pf, as well as the number of scheduled users per cell K, that 
maximize the overall spectral efficiency within the cell. In the 
following, we obtain expressions for the spectral efficiency 
adopting MRC and P-ZFC combining. Without loss of gener¬ 
ality, we will hereafter assume that the user coordinates zik 
are uniformly distributed within the cells. 


B. Achievable Spectral Efficiency 


We first give the following preliminary dehnitions. Let I; 
be the set of the PfK users in cell I G C with the highest 
d/(-); likewise, let £i be the set of the — Pf)K users in cell 
I G C with the lowest dp-). Furthermore, we introduce the 
functions 





( 6 ) 

(7) 


for 7 = 1,2, where the ratio dj{zik)/dpzik) expresses the 
relative strength of the interference received at BS j from user 
k in cell 1. Lastly, let C, c C denote the set of cells using the 
same subset of pilots as £j^ 


Theorem 1. The ergodic achievable spectral efficiency in cell 
j, when the users are uniformly distributed within the cell, is 
given by 

SE,4b(i-|) (8) 

X (/3f log2(l + SINRxJ + (1 - ^f)log2(l -f SINRxj) 

with B defined in © and where the expressions of the signal- 
to-interference-plus-noise ratios (SINRs) are given in (0-([Tg 
for MRC and in ( |1 l| i-( [T^ for P-ZFC, respectively, at the top 
of the page. 


'^Observe that Cj = C occurs only when /? = 1. 


Proof: The expression for the ergodic achievable spectral 
efficiency in cell j without FPR is given in pT| Th. 1]: 

SE, 4 (l-_)^log2(l + SINR,fe) (13) 

k=l 

where SINR^x is the SINR corresponding to user k in cell j 
derived in Th. 1, 2] for MRC and P-ZFC and 



We first note that, if the coordinates of all users within a 

('y) 

cell have the same distribution, then p - f loses its (implicit) 
dependence on the user index k and, as a consequence, we 
can rewrite the following term as 

E E = E/^5? E ^ E ^5? 

l£C m—l l£C m—l iGCj 

(15) 


where the last equality derives from the fact that 
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= 




B, 

0 , 


if I G Cl¬ 


ot her wise. 

By extending (|T5]l to the case of FPR, we can write 
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Figure 3. Comparison between pilot reuse (PR) |1Q| and fractional pilot reuse 
(FPR) with MRC: (a) spectral efficiency of the cell as a function of TV; (b) 
Optimal number of scheduled users as a function of TV. 
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Figure 4. Comparison between pilot reuse (PR) |10| and fractional pilot reuse 
(FPR) with P-ZFC: (a) spectral efficiency of the cell as a function of TV; (b) 
Optimal number of scheduled users as a function of TV. 


and we obtain the SINRs in Finally, exploiting the 

fact that \Ij\ = f3{K, we arrive at the expression of the spectral 
efficiency in ([^. ■ 

It is straightforward to conclude that, when (3{ = 0, we have 
1/ = 0 and the spectral efficiency in •S' collapses to the 
expression given in CD. 

Remark 1. The expressions above can be easily extended 
to the case where the K mobile users are not uniformly 
distributed over the cell at the cost of less tractability. More 
specifically, when the average relative strength of a user 
depends on its index, then the terms and cannot be 
removed from the inner summation. Nonetheless, the resulting 
formulas can be evaluated by Monte-Carlo simulations if 
required. 

Now, let us study the performance of the proposed scheme 
based on FPR when N ^ oo. 


Corollary 1. The spectral efficiency in in the large-antenna 
regime (N —> oo), when the users are uniformly distributed 
within the cell, becomes 


= /f(i - f) (i+ 


1 + 


iec\{3} 
1 




( 21 ) 




with B defined in 0. 

Evidently, the asymptotic spectral efficiency derived in Corol¬ 
lary does not depend on the type of combining employed 
at the BS. Again, we observe that, when /3f = 0, the same 
expression of the asymptotic spectral efficiency given in CD 
Cor. 3] is obtained. 

Unfortunately, as in the expressions obtained in ([^ and 
in ( |2 T] i do not allow for closed-form extraction of the optimal 
{/3,/3f,iC}; moreover, a further complication of our scheme 
is that the terms ^nd intrinsically depend on /3f. 

Therefore, in the next section, we optimize the system param¬ 
eters using numerical evaluation for a particular scenario. 

IV. Numerical results 

In this section, we assess the performance of the proposed 
scheme based on FPR through numerical evaluation. In doing 
so, we adopt an equivalent setup to 03 and compare our 
results against those obtained without using FPR. 

Let us consider a symmetric, infinitely large network of 
hexagonal cells with radius r > 0. Furthermore, let us assume 
a standard pathloss model in which the variance of the channel 
attenuation is given by dj(z) = ||z — bj||^”, where k > 2 is 
the pathloss exponent and bj G specifies the coordinates 
of BS j: in this context, we have (cf. @-0) 

djjzik) _ / \\zik - b ;||2 
di{zik) V||z;fc - bj|| 2 / 


-f (1 - /3f) log 2 


( 22 ) 
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Figure 5. Spectral efficiency with N = 1000 antennas as a function of /3f. 




MRC 

P-ZFC 

N = 

10® 

5.6% 

93.1% 

N = 

10^ 

10.4% 

125.8% 

N = 

10® 

21.2% 

117.4% 

N = 

10^ 

28.6% 

81.3% 


Table I 

Gains of FPR with respect to (To). 


For our numerical simulations, the same parameters as in 
m are used: we take into account the interference produced 
by 3 tiers of surrounding cells, a coherence block of T = 1000 
channel uses, a pathloss exponent k = 3.5, an average SNR 
between any user and any antenna at the corresponding BS 
of pjcP' = 10 dB, and a number of BS antennas in the range 
N S [10,10^]. Lastly, the terms and are computed 
by Monte-Carlo simulations as follows. For each value of K, 
we generate 10® instances of K uniformly distributed user 
coordinates in every cell I with a minimum distance of 0.14r 
from the corresponding BS; then, for each value of /3f, i) the 
PfK coordinates that are closest to the BS are assigned to I; 
and is computed as in and ii) the other coordinates 
are assigned to £i and is computed as in Q. 

Remarkably, using FPR always improves the spectral ef¬ 
ficiency with respect to fTO) , as summarized in Table 
Figure [^a) and Figure [3[b) show the spectral efficiency 
and optimal number of scheduled users, respectively, when 
using MRC combining; despite obtaining quite similar optimal 
values of K with respect to the baseline flO) , we observe that 
adopting FPR allows to anticipate the switching point from 
/3 = 3 to /3 = 1, i.e., the channels become quasi-orthogonal 
at lower values of N. Likewise, Figure illustrates the same 
results for the case of P-ZFC, which produces significantly 
improved spectral efficiencies with respect to MRC. Note that, 
with P-ZFC, choosing /3 = 1 yields a higher spectral efficiency 
with respect to /3 = 3 even for low values of N, which 
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Figure 6. Spectral efficiency using P-ZFC and /3 = 1 as a function of /3f for 
different values of TV. 

always allows to schedule more users than without FPR (see 
Figure Qb)). 

Let us more generally compare MRC and P-ZFC: in the 
former, the impact of our approach is more significant as the 
number of antennas at the BS increases whereas, for the latter, 
our method exceptionally doubles the spectral efficiency in 
the whole range of N under exam. In addition, by comparing 
Figure [^a) and Figure [^a), it emerges that FPR with MRC 
always outperforms the baseline ng with P-ZFC. 

As expected, the spectral efficiency in ([^ is concave in /3f 
for any fixed K. The spectral efficiency as a function of is 
depicted in Figure]^ with N = 10® BS antennas and using the 
optimal value of K for each setup: in particular, we observe 
that the optimal /3f changes with both the integer reuse factor 
/3 and the combining technique. Lastly, the spectral efficiency 
obtained with P-ZFC /3 = 1, i.e., the best combination of 
integer reuse factor and combining technique, is illustrated as 
a function of /3f in Figure for some values of N. 

V. Conclusions 

In this paper, we propose a novel pilot allocation scheme 
for multi-cell massive MIMO systems, where users close to 
their respective base stations get to reuse the same set of pilot 
sequences across the whole system. We derive expressions 
for an achievable rate in this scenario and obtain both the 
pilot reuse parameters and the number of scheduled users per 
cell that maximize it. Results show a significant gain in terms 
of spectral efficiency with respect to the existing pilot reuse 
baseline. 
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